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A new measure to identify a small-scale dissipation region in collisionless magnetic reconnection 
is proposed. The energy transfer from the electromagnetic field to plasmas in the electron's rest 
frame is formulated as a Lorentz-invariant scalar quantity. The measure is tested by two-dimensional 
particle-in-cell simulations in typical configurations: symmetric and asymmetric reconnection, with 
and without the guide field. The innermost region surrounding the reconnection site is accurately 
located in all cases. We further discuss implications for nonideal MHD dissipation. 

PACS numbers: 52.35.Vd, 94.30.cp, 95.30.Qd, 52.27.Ny 



Magnetic reconnection [1 is a fundamental process in 
many plasma systems, ranging from laboratory and solar- 
terrestrial environments to extreme astrophysical set- 
tings. The violation of the ideal condition, E^v xB ^ 0, 
is essential to allow the magnetic flux transport across the 
reconnection point. The critical "diffusion region" (DR) 
where the ideal condition is violated is of strong interest 
for understanding the key mechanism of reconnection. In 
collisionless plasmas, since ions decouple first from the 
magnetic fields, it is thought that the DR consists of an 
ion-scale outer region and an electron-scale inner region. 

In two-dimensional (2D) reconnection problems in the 
x-z plane, a popular criterion to identify the innermost 
"electron diffusion region" (EDR) is the out-of-plane 
component of the electron nonideal condition, ^ 0, 
where 
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and the electron pressure tensor. In particular, it is 
known that the divergence of the pressure tensor sustains 
a finite Ey = Ey at the reconnection point, arising from 
local electron dynamics [H [3] . 

Recent large-scale particle-in-cell (PIC) simulations 
have shed light on the electron-scale structures around 
the reconnection site. Earlier investigations [H [5] found 
that the EDR identified by 7^ [4 or the out-of-plane 
electron velocity [5 extends toward the outflow direc- 
tions. Previous research has suggested that the EDR has 
a two-scale substructure: the inner EDR of E'* > and 
the outer EDR of < with a super- Alfvenic electron 
jet in [7]. Satellite observations found similar signatures 
far downstream of the reconnection site [8^. The roles of 
these EDRs are still under debate; however, there is a 
growing consensus that only the inner EDR or a similar 
small-scale region should control the reconnection rate 
[71|9j[T0j. Importantly, it was recently argued that the 
outer EDR is non- or only weakly dissipative, because the 
super- Alfvenic jet and Ey ^ condition stem from pro- 
jections of the diamagnetic electron current in a suitably 
rotated frame [10]. 

Meanwhile, a serious question has been raised by 



numerical investigations on asymmetric reconnection, 
whose two inflow regions have different properties such 
as in reconnection at the magnetopause [11] |T2] . It was 
found that various quantities including E* fail to locate 
the reconnection site in asymmetric reconnection, espe- 
cially in the presence of an out-of-plane guide field [TT] . 
Considering the debate on inner or outer EDRs and the 
puzzling results in asymmetric reconnection, it does not 
seem that £^*^0 is a good identifier of the critical region. 

In this Letter, we propose a new measure to identify 
a small, physically significant region surrounding the re- 
connection site. We construct our measure based on the 
following three theoretical requirements. First, we are 
guided by the notion that dissipation should be related 
to nonideal energy conversion. Second, we desire a scalar 
quantity. If we use a specific component of a vector, we 
have to choose an appropriately rotated frame [10 . Using 
a scalar quantity instead, we do not need to find the right 
rotation in a complicated magnetic geometry. Third, it 
should be insensitive to the relative motion between the 
observer and the reconnection site. For example, the re- 
connection site can retreat away [13 , or, for example, the 
entire reconnection system may flap over a satellite due 
to the magnetospheric motion. 

Our strategy is as follows. We choose a frame that 
can be uniquely specified by the observer. Among several 
candidates, we choose the rest frame of an electron's bulk 
motion because it would be the best one to characterize 
electron-scale structures. Next we consider the energy 
transfer from the field to plasmas in this frame, which is 
a scalar quantity. We then expand it with observer-frame 
quantities. The obtained measure meets all three re- 
quirements. It is a Lorentz invariant (frame-independent 
scalar) and is related to the nonideal energy transfer. 

We follow the spacelike convention (—,+,+,+). Let us 
start from the electromagnetic tensor F^^, 
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Using a 4- velocity {u^) = 7(0, v), where 7 is the Lorentz 
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factor 7 = [1 — {v/c)'^] we obtain a 4- vector of the 
rest-frame electric field [14^ , 

gM = AMg/^, (3) 

Here the prime sign ' denotes the properties in the rest 
frame of an arbitrary motion ^ and A is the inverse 
Lorentz transformation from the moving frame. The 
components of and e'^ are given by, 

(eA') = (2^, + vxB)), (e'^) = (0, E'). (4) 

We also use the 4-current (J^) = {pcC^j)^ where pc is the 
charge density. The current can be split into the conduc- 
tion current and the convection current, a projection 
of the motion of the non-neutral frame [15 , 

= jf^ + p/^^/i = jf^ + c-\-ru^)u^, (5) 

such that (/^) = (0, j') is purely spacelike. 

Let us define a dissipation measure D, the energy con- 
version rate in the moving frame. The contraction of the 
covariant and contravariant vectors gives us a Lorentz- 
invariant scalar, 

= j[j-{E^vxB)-p,{vE)]. (6) 

Choosing the frame of electron bulk motion (the number 
density's flow), we obtain the electron-frame dissipation 
measure^ 

De = 7e [j -{E^VeXB)- p,{Ve ' E)] . (7) 

In the nonrelativistic limit, one can simplify Eq. ^ by 
setting 7e ^ 1. One can confirm this by multiplying 
j' = qniv[ = (j - pcVe) and E' = E* = {E ^VeX B). 

In ion-electron plasmas, since = Qi^iU^ — ng'U^), 
where n is the proper density, we obtain the following 
relation between the electron-frame and ion-frame mea- 
sures, 

rieDe = UiDi. (8) 

Such a symmetric relation is reasonable, as ions are the 
current carrier in the electron's frame and vice versa. 
If ions consist of multiple species, UeDe = ZgUsDs^ 
where s denotes ion species and Z is the charge number. 

To see how our measure characterizes the reconnection 
region, we have carried out 2D nonrelativistic PIC simu- 
lations. The length, time, and velocity are normalized by 
the ion inertial length di = c/ujpi^ the ion cyclotron fre- 
quency and the ion Alfven speed c^i, respectively. 
The mass ratio is m^/me = 25, and the electron- ion tem- 
perature ratio is T^/Ti = 0.2. Periodic (x) and conduc- 
tive wall (z) boundaries are used. Four runs (1-4) are car- 
ried out. Runs 1 and 2 employ a Harris-like configuration, 
B{z) = Bot8inh{2z)x and n{z) = no[0.2 + cosh"^(2^)]. 



The domain of [0,102.4] x [-25.6,25.6] is resolved by 
1600^ cells. 2.6x10^ particles are used. The speed 
of light is c = 10. In run 2, we impose a uniform 
guide-field By = Bq. Runs 3 and 4 employ asym- 
metric configuration. Since no kinetic equilibrium is 
known, we employ the following fiuid equilibrium pro- 
posed by Ref. [16 , B{z) = Bo[^ + tanh(2z)] x and 
n{z) = no[l — I tanh(2z) — | tanh^(22;)] . Across the 
current sheet, magnetic fields and the density vary from 
—Bo/2 and no to 3Bo/2 and no/3. The domain of 
[0, 64] X [-12.8, 12.8] is resolved by 1000x800 grid points. 
9x10^ particles are used. The speed of light is c = 20. 
In run 4, a guide-field By = Bq is added. In all runs, 
reconnection is triggered by a small fiux perturbation. 
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FIG. 1. (color online). Snapshots of run 1 at t = 60, aver- 
aged over . (a) The nonideal electric field Ey and (b) the 
electron- frame dissipation measure De (Eq. 0). 



The panels in Fig. [T] present the popular measure Ey 
and the electron-frame dissipation in run 1 in the 
well-developed stage. They are normalized by caiBq and 
CAiBojoi respectively. Another option is to employ the 
upstream normalization [6l |7] or its hybrid extension for 
asymmetric cases [T7|, but these are beyond the scope 
of this paper. All quantities are averaged over to 
remove noise. In Fig. [ija), one can recognize a positive 
Ey region near the reconnection site and a negative Ey 
channel which extends to the outfiow direction. They 
correspond to the inner and outer EDRs 1^ . On the 
other hand. Fig. gives a different picture. There is a 
positive De region near the reconnection site, indicating 
that the strong energy transfer occurs there. Hereafter 
we call it "dissipation region" of I^e > 0. At the recon- 
nection point, a main contributor to j' ■ E' is jyEy. In 
this case, the charge term pcVe-E is responsible for —25% 
of the total value, due to significant charge separation. 
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l^i — ^e|/(^i +^e) ^ 15%. As wc move to the outflow di- 
rection, jyEy is gradually replaced by jxEi^- This makes 
the dissipation region longer than the inner EDR. Based 
on the scale height, the aspect ratio of the dissipation 
region, 9.5 : 0.72, is similar to the universal reconnection 
rate of 0.1. Outside of minor fluctuations, there are no 
signiflcant structures in the downstream region {x > 55). 

We have also studied the ion-frame dissipation Di in 
this case. The spatial proflle closely resembles to that 
of as indicated by Eq. ([8| for such a quasineutral 
plasma (ng — rii). 
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FIG. 2. (color online). The dissipation measure De (Eq. 0) 
at t=60 in other runs: (a) run 2, (b) run 3, and (c) run 4. 
The dash line show the field reversal line, Bx — 0. 

The panels in Fig. [2] show the dissipation measure 
in three other runs. One can see that it excellently iden- 
tifles compact regions about the reconnection sites in all 
cases. In the symmetric guide-fleld run [Fig. 2(a)], the 
dissipation region is tilted slightly anticlockwise. This 
is associated with the electron cavities with a parallel 
electric fleld along one pairs of separatrices [18 . In the 
asymmetric cases, De excellently works even in the most 
challenging case with a guide-fleld [11 . The fleld reversal 
lines {Bx = 0), shown by the dashed lines, are located 
inside our dissipation regions. The peak amplitude of 
De is high in the guide fleld cases, as the strong electron 
current is conflned. This deserves further investigation, 
because the kinetic dissipation mechanism is different in 



guide-fleld cases [3^. There are charge-separated regions 
along the separatrices in the guide-fleld cases and on the 
boundary with the upper inflow regions in the asymmet- 
ric cases. Such charge-separation effects are included in 
De via the last term in Eq. Q, which usually improves 
the identiflcation of the dissipation region. 
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FIG. 3. (color online). Composition of the reconnection elec- 
tric field Ey (Eq.[l]) along the infiow line at (a) x = 48 in run 
1 and (b) x = 31.5 in run 3. The shadow presents a rescaled 
value of De to indicate the dissipation region. 

Let us focus on the dissipation region in runs 1 and 
3. The panels in Fig. [3] show the composition of the 
reconnection electric fleld Ey along the inflow lines, based 
on Eq. ([T]). The shadows indicate the rescaled amplitude 
of the dissipation measure De. 

In the symmetric case [run 1; Fig. 3(a)], Ey is bal- 
anced by the electron pressure tensor term [2] at the cen- 
ter and the bulk inertial term in surrounding regions. 
The dissipation region is located in a narrow region on 
local electron inertial scale, de^O.5-6. There, the bounc- 
ing electrons carry a strong electron current, and jy is 
intense accordingly. We also note that the electron ideal 
condition is weakly violated on the larger scale of the lo- 
cal ion inertial length, ^2-3. Such an outer structure is 
related to ion's decoupled motion, and we will see a clear 
two-scale structure at sufficiently large times [19 . 

In asymmetric reconnection, it is known that the 
fleld reversal {Bx = 0) and the flow stagnation points 
{vizi^ez = 0) are usually not collocated [I6l [TTl [20] . In 
our run [Fig. 3(b)], the electron flow stagnation point 
is located on the upper side {z :^ 0.63) while the fleld 
reversal is on the lower side {z —0.3). Therefore 
the motional electric fleld even becomes negative around 
the dissipation region. Importantly, Ey remains positive 
there and it resonates with a main current jy > 0. In- 
terestingly, jzE'^ also contributes to De in the lower side 
{z ~ —0.3). On the other hand, even though E'^ is an 
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order of magnitude larger than Ei^ and Ey on the up- 
per side {0 < z < 0.5) [iT, it does not contribute to 
because the vertical current jz is negligible. Physically, 
such a strong E'^ is overemphasized by the diamagnetic 
drift where the density gradient is strong. Regarding Ey^ 
we find that the ideal condition is violated in the lower- 
side upstream of z < 0. The pressure tensor terms, both 
dxPexy and dzPeyz^ are not negligible there. We find that 
they stem from the gyrotropic electron pressure tensor 
^ e ^ Pe±I-^{Pe\\—Pe±)BB/B'^ in the upstream region. 
This tells us that these effects are due to the drift mo- 
tion of gyrating electrons. We also notice that the time 
derivative term (dt) is a key contributor on the upper side 
{z ~ 0.6), because the reconnection site moves upward 
very slowly. Even when the structure is stationary in a 
frame d/dt = 0, the time derivative dt = — (t^frame • V) is 
not always zero in the observer frame. The fact dt 
renders the analysis more difficult and it provides another 
motivation for using a frame-independent measure. 

Let us discuss the relevance for the MHD energy bud- 
get. For simplicity we limit our discussion to the non- 
relativist ic regime. Defining an MHD velocity i^mhd = 
{miTiiVi + merieVe) / {miTii + mgne), we find 
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The total energy transfer can be decomposed to 



j ■ E = {j X B) - Vmhd + PcE • Vmhd + D 
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The first term stands for the work done by the Lorentz 
force on the MHD fiuids. This operates in the ideal MHD 
also. The second term is the work done by the Coulomb 
force, also interpreted as the energy transfer by the con- 
vection current. These two disappear in the MHD frame. 
The last term is responsible for the nonideal energy trans- 
fer. In a quasineutral plasma rii ^ the positive I^mhd 
(c± De) plays the same role as an irreversible dissipation 
given by rjj^ > in the resistive MHD. One can see that 
positive De regions are enhanced and localized around 
the reconnection sites in Figs, [ijb) and [5] 

From the kinetic viewpoint, gyrating particles undergo 
various drift motions such as V5, diamagnetic, and cur- 
vature drifts. When non-£^ x B drifts appear, the bulk 
fiow no longer comoves with the field lines, and there- 
fore the ideal condition is no longer a useful concept. 
Drift motions lead to the electromagnetic energy dissi- 
pation, if and only if they involve nonideal energy con- 
version. For example, we do not recognize a significant 
nonideal energy transfer in the outer EDR in Fig. [ijb), 
where E^y^O due to the diamagnetic effect [lOj. On the 
other hand, around the reconnection sites, nongyrotropic 
or field-aligned electrons carry intense currents and then 
they enhance the nonideal energy transfer. 

In summary, we have proposed an electron-frame dis- 
sipation measure De [Eq. ^] to identify a physically rel- 
evant, small-scale region surrounding the reconnection 



point. We have demonstrated that it works excellently 
in typical configurations. Furthermore, we identified its 
relation to nonideal MHD dissipation, which is essential 
to the reconnection problem. 

Our finding will benefit several research fields. One 
is the satellite observation of reconnection. NASA is 
preparing the Magnetospheric Multiscale (MMS) mission 
to observe the electron-scale structures in near-Earth re- 
connection sites. By using Eq. ^ with 7e ^ 1, one 
can identify the dissipation region regardless of the mo- 
tion and the orientation of the reconnection site. An- 
other example is relativist ic astrophysics. Owning to a 
growing attention to reconnection, numerical modeling of 
relativistic reconnection has been growing in importance 
[2T] , but basic properties are much less known than in 
the nonrelativistic case. Our Lorentz- invariant measure 
can be readily applied to the relativistic dissipation re- 
gion problem [22 . Further numerical work is desirable 
to further test our measure in three dimensions for these 
systems. 
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